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The nature of solar wind turbulence in the dissipation range at scales much smaller than the large MHD scales 
remains under debate. Here a two-dimensional model based on the hybrid code abbreviated as A.I.K.E.F. 
is presented, which treats massive ions as particles obeying the kinetic Vlasov equation and massless elec- 
trons as a neutralizing fluid. Up to a certain wavenumber in the MHD regime, the numerical system is 
initialized by assuming a superposition of isotropic Alfven waves with amplitudes that follow the empirically 
confirmed spectral law of Kolmogorov. Then turbulence develops and energy cascades into the dispersive 
spectral range, where also dissipative effects occur. Under typical solar wind conditions, weak turbulence 
develops as a superposition of normal modes in the kinetic regime. Spectral analysis in the direction parallel 
to the background magnetic field reveals a cascade of left-handed Alfven/ion-cyclotron waves up to wave 
vectors where their resonant absorption sets in, as well as a continuing cascade of right-handed fast-mode and 
whistler waves. Perpendicular to the background field, a broad turbulent spectrum is found to be built up of 
fluctuations having a strong compressive component. Ion-Bernstein waves seem to be possible normal modes 
in this propagation direction for lower driving amplitudes. Also signatures of short-scale pressure-balanced 
structures (very oblique slow-mode waves) are found. 



PACS numbers: 94.05.Lk, 94.05.Pt 

I. INTRODUCTION 

The solar wind is a dilute plasma and known to be in 
a highly turbulent state 1 2 . It exhibits fluctuations in the 
electromagnetic field, the plasma density, and bulk flow 
velocity over a wide range of scaled. However, the nature 
of solar wind turbulence in the intermediate wavenumber 
regime, situated between the large inertial MHD scales 
and the small dissipative electron scales, is not well un- 
derstood. Especially, the role of oblique wave propaga- 
tion with respect to the background Parker field is cur- 
rently under debate. Some authors favor a more or less 
independent behavior of the so called slab component 
(parallel with respect to the background magnetic field) 
and 2D turbulence (perpendicular) 4 6 . While in this pic- 
ture the highly oblique 2D-turbulence is believed to be 
an example of a strongly turbulent plasma state with 
high-order correlations between the fluctuating quanti- 
ties, the slab component is assumed to be describable 
within the framework of weak turbulence theory, i.e. as a 
superposition of normal modes. Other authors interpret 
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the observed anisotropy as being the result of a kinetic 
Alfven wave (KAW) cascade^, which would suggest a 
preference of fluctuations with k± ^> kn . 

Incompressible, dissipative MHD simulations in two 
dimensions 9 and three dimensions^ show the develop- 
ment of an anisotropy in turbulent spectra with respect 
to the background magnetic field. In these cases, an ini- 
tially isotropic spectrum is shown to evolve towards an 
anisotropic configuration with preferred wavevectors per- 
pendicular to the background field. These authors have 
shown that the strength of dissipation in the simulations 
is crucial for the development of anisotropy. The theoret- 
ical basis of the geometry in the presence of a strong back- 
ground field was discussed by Montgomery and Turner 4 
in the framework of a perturbation analysis of the in- 
compressible MHD equations. The resonant conditions 
ki + &2 = ks for the wavevectors and uj\ + 002 = ^3 for 
the wave frequencies in a resonant three- wave interaction 
were found to favor a perpendicular cascade and lead to 
a suppression of energy transfer in the parallel direction. 
The strongest turbulent coupling is found to occur if all 
three frequencies are about zero, a situation which is clos- 
est to the hydrodynamic case. This can be completely 
achieved only if the three-dimensional nature of the solar 
wind fluctuations is fully resolved, like with a 3-D code 
that permits to consider both dimensions perpendicular 
to the background field. However, it was claimed that 
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higher-order couplings could in fact also yield a parallel 
cascade, yet with a slower growth. These higher-order 
couplings, the dispersion of the normal modes, and pos- 
sible compressive effects lead to additional complications 
in the kinetic regime and make a numerical treatment 
necessary. 

Measurements made by the four Cluster spacecraft 
have provided further insights into the nature of so- 
lar wind turbulence, since multi-spacecraft detections of 
magnetic fluctuations even permit the analysis of their 
three-dimensional dispersion properties. However, these 
measurements support different interpretations. A fully 
evolved nonlinear turbulent state with a preferred 2D- 
component beyond the MHD range was found by Alexan- 
drova et al. 11 The dispersion relation of this mainly per- 
pendicularly structured component might reflect a non- 
linear cascade that also occurs at small scales. However, 
an interpretation on the basis of normal modes, such 
as the right-handed circularly polarized fast /whistler 
(F/W) waves, is still possible^. 

There is strong evidence that the field-parallel com- 
ponent consists at least partly of Alfven/ion-cyclotron 
(A/IC) waves, which are dispersive left-hand circularly 
polarized electromagnetic normal modes of a plasma. 
The temperature anisotropics and beam structures ob- 
served in the solar wind proton distribution functions 
were explained as resulting from cyclotron-resonant inter- 
actions of the ions with these waves^HUl Recently, direct 
wave measurements have confirmed the existence of A/IC 
waves in the solar wind 16 . Yet, due to cyclotron-resonant 
wave-particle interactions, they are strongly damped at 
wavenumbers corresponding to the inverse inertial length 
of the resonant ions^, and thus they are not expected 
to exist at higher wavenumbers. A coexistence of left- 
handed and right-handed modes has been recently sup- 
ported by measurements of the angle distribution of the 
magnetic helicity in the solar wincP^. Possible candi- 
dates for normal modes beyond the resonant wavenum- 
ber range are the F/W modes, which may remain af- 
ter the dissipation of the A/IC waves^Hl, or an ongoing 
cascade^ of dispersive KAWs. Both these wave modes 
are right-hand polarized under the conditions prevailing 
in the solar wind. It is observed that the right-hand po- 
larized waves survive the spectral break, which indicates 
the transition from the inertial range to the ion dissipa- 
tive scales, and that they can exist at higher wavenum- 
bers without damping, until the resonant electron scales 
are finally reached 20 . 

There are other numerical hybrid models available, 
which treat the kinetic cascade of turbulence beyond the 
MHD scales. For example Parashar et alP^ achieve an 
efficient turbulence driving by assuming an Orszag-Tang 
vortex in the plane perpendicular to the background mag- 
netic field. This configuration leads to perpendicular ion 
heating, which seems to be not connected to cyclotron- 
resonant effects. In a later workP^l, the authors also in- 
cluded electron pressure effects and found indications for 
Bernstein waves. However, the spectrum did not show 



significant presence of Alfven or whistler waves but rather 
a dominance of zero-frequency structures for their geom- 
etry and in the strongly turbulent regime. These simula- 
tions show that the presence of the electron pressure term 
in the generalized Ohm's law can substantially modify 
the properties of the modeled turbulence. A similar sit- 
uation has recently been treated in another worW 2 ^ with 
a special focus on the dissipation of the energy. The for- 
mation of intermittent current sheets on short time scales 
and the related demagnetization of the plasma ions lead 
to a stochastic heating 2 -^. 

Our numerical simulation work is focused on study- 
ing the transition of isotropic MHD turbulence to non- 
isotropic kinetic fluctuations on intermediate scales. For 
this purpose, at least two-dimensional numerical simula- 
tions are necessary, by which one can analyze the evo- 
lution of turbulence in different directions with respect 
to the constant background magnetic field. The model 
system is initialized by a superposition of linear MHD 
waves. No further external or ongoing driving force is 
applied, and thus the system will evolve freely from its 
initial state. Two separate simulation runs are presented, 
which only differ in the amplitude of the initial pertur- 
bation. 



II. NUMERICAL METHOD 

The A.I.K.E.F. code is a numerical hybrid code, which 
treats ions as particles following the characteristics of 
the Vlasov equation and electrons as a massless charge- 
neutralizing fluid. It is based on the work by Bagdonat 
and MotschmanrP^, and its basic properties have later 
been described extensively by Bagdonat 26 . The code was 
completely revised, and the possibility for adaptive mesh 
refinement was included by Miiller et alP^The equations 
of motion for a single proton are 



for the velocity v p and the spatial location x p . The force 
acting on the particle with charge q p and mass ra p is the 
Lorentz force, which is due to the electric field E and 
the magnetic field B. The speed of light is denoted by c. 
The electron fluid equation of motion delivers the electric 
field as 

^ = --ix5-— Vp e , (3) 
c n e e 

where the electron bulk velocity is denoted by u ei the 
electron number density by n e , the elemental charge by 
e, and the electron fluid pressure by p e . The electrons are 
assumed to be isothermal, and thus the pressure depends 
on the electron number density according to p e oc n e . 
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Quasi-neutrality requires that n p = n e . The proton den- 
sity and the proton bulk velocity u v are obtained as the 
first two moments of the proton distribution function at 
every step in time. 

The magnetic field is obtained from the induction 
equation following from Faraday's law: 



dB 

~dt 



= V x (u p x B^j - V x (|^-V x B x B^j . 



(4) 



The ion charge density g c corresponds to the proton num- 
ber density calculated as the zeroth velocity moment of 
the distribution function of all protons. 

The boundaries of the simulation box are set to be pe- 
riodic, and the particles are initialized with a Maxwellian 
velocity distribution that is shifted to the given val- 
ues for the initial bulk velocities. The width of the 
Maxwellian distribution is determined by the respective 
species' plasma beta, which represents the ratio of the ion 
thermal to the magnetic energy density. The beta is set 
to /3 p = 0.05 for the protons. The electron beta is fixed 
at /3 e = 0.5. In this beta regime, even low amplitudes 
of magnetic fluctuations will have a strong influence on 
the motion of the particles due to their high magnetiza- 
tion. All spatial length scales are normalized and given 
in units of the proton inertial length, i v = c/a; p , with the 

proton plasma frequency cu p — ^4irn p q^/m p . All time 

scales are in units of the inverse proton gyro-frequency 
Q p = q p B/(m p c). In these units, the two-dimensional in- 
tegration box has a size of 250 x 250, which is covered by 
1024 x 1024 cells, each of which filled with 500 superparti- 
cles representing the real number density of the protons. 
In this spatial configuration, one grid step corresponds 
to about 0.24^ p in the calculation. The data is recorded 
with a grid resolution of Ax = 0.5£ p to keep the amount 
of data reasonable. This corresponds to a maximum re- 
solvable wavenumber of £ p k res « 7r/0.5 « 6.3, which is 
adequate for an appropriate coverage of the transition 
from MHD into the kinetic regime. For all vector quan- 
tities, three components are evaluated. Approaches of 
this form are sometimes referred to as 2.5D simulations. 
A divergence-cleaning algorithm is applied to guarantee 
numerical stability. 

The initial magnetic field is given as a superposition of 
linear Alfven waves according to 
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SB X = ^2 ^2 frnCos(/c n ?/sin?9 m + k n z cos <d m + p n>m ), 

m=l n=l 

(5) 

where the constant background field Bq is aligned along 
the z-axis. A random phase shift p n ,m is applied to each 
wave. The amplitudes b n are fixed in such a way that 
the power spectrum follows a Kolmogorov power law^ in 
wavenumber with the scaling oc /c -5 / 3 . The total power 
of the composed wave field is made equal to the power of 
a monochromatic wave with SB = 0.01£?o for Run A and 
to 5B = O.l^o for Run B. The first value corresponds to a 



situation with a weaker turbulence level or a higher back- 
ground field as prevailing in the solar corona, while the 
second value represents realistic solar wind conditions. 
The angles # m cover 360° of propagation directions by 
m max = 50 discrete values, and the spectrum ranging be- 
tween fco =0.05 and /c max = 0.2 is covered by n max = 20 
waves. The upper limit & max is quite high compared to 
typical MHD scales but, still, in the dispersion-less range 
in first order. The initial velocity is obtained from the 
Alfvenic polarization relation 
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with the Alfven speed Va = £ p /n p . 



III. RESULTS 

In this section some results of the numerical simulation 
runs are shown. We ran the code for a sufficiently long 
time so that an evolved nonlinear dynamic plasma state 
could be expected. This was typically the case after an 
evolution time of about 500 gyro-periods, at which time 
the system was analyzed. 



A. Results for simulation Run A 

For the analysis of the results from Run A, Fourier 
transformations in two dimensions were applied subse- 
quently to the magnetic field data and the density data. 
The power spectral density was then calculated, and 
could be shown to depend on the wavenumbers k z for the 
direction parallel to the background field and k y perpen- 
dicular to it. The resulting magnetic field power spectral 
density is shown in Fig. [I] and the power spectral density 
of the compressive fluctuations is shown in Fig. [2] 

Apparently, the magnetic field fluctuations show a pre- 
ferred alignment with the direction perpendicular to the 
background magnetic field. Especially, turbulence energy 
is spreading in this direction to much higher wavenum- 
bers than in the parallel direction. But also the parallel 
fluctuations at wavenumbers beyond the initialized range 
are excited and gain energy. The compressive density 
fluctuations, however, are mainly aligned perpendicular 
to the background field and have almost no components 
parallel to Bo. 

Cuts through the two-dimensional spectra along the 
perpendicular direction are shown in Fig. [3] The power 
spectral density of the compressive fluctuations in the 
proton number density is enhanced for k y > 1 and fol- 
lows mainly a power law with a slightly steeper index 
than /c -5 / 3 in this range. The magnetic field spectrum is 
flatter in the dispersive range. 

To study the nature of these fluctuations Fourier trans- 
formation can also be used, but here is applied in the 
time domain leading to the corresponding dispersion 
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FIG. 1. Two-dimensional power spectral density of magnetic 
field fluctuations in arbitrary units for Run A. The back- 
ground magnetic field is oriented along the /c z -axis. The cas- 
cade of energy to higher wavenumbers occurs preferentially in 
perpendicular direction. 



FIG. 4. Dispersion relation of magnetic field fluctuations par- 
allel to the background magnetic field at t ~ 500 for Run A. 
The overplotted red lines indicate the F/W (upper line) and 
A/IC (lower line) wave branches calculated directly from the 
cold-plasma dispersion relation. 
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FIG. 2. Power spectral density of compressive fluctuations in 
the proton number density in arbitrary units for Run A. The 
spatial variation of the fluctuations occurs mainly perpendic- 
ular to the background magnetic field. 
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FIG. 3. Power spectral density of fluctuations in the magnetic 
field and density at t — 500 along the perpendicular direction 
for Run A. A power law with spectral index —5/3 is shown 
additionally, in order to estimate the slope of the spectrum. 



diagrams. Therefore, first the two-dimensional spatial 
Fourier transformation is applied, and the result is taken 
in one dimension only (parallel or perpendicular to Bo), 
and then the data are again Fourier transformed yet in 
time. The result for the parallel magnetic field dispersion 
is shown in Fig. [4j 

Most of the power is found at low wavenumbers and 
low frequencies and still close to the initial distribution of 
the waves. Spreading of power at low parallel wavenum- 
bers with low frequencies is an indication for ongoing 
MHD turbulence. Two sharp branches can clearly be 
seen beyond the initial wavenumber limit at /c max = 0.2. 
For an easier identification of them, the cold-plasma dis- 
persions for the left-handed A/IC waves and for the right- 
handed F/W waves are additionally shown in the paral- 
lel dispersion diagram 29 . The observed parallel disper- 
sion agrees well with that of the linear normal modes. 
The A/IC branch ends at a wavenumber value below 
1, whereas the right-handed branch continues to higher 
wavenumbers and frequencies. 

The perpendicular magnetic field dispersion is depicted 
in Fig. [5] Since the density fluctuations are mainly per- 
pendicularly oriented, their parallel dispersion diagram 
is not shown here. But the perpendicular dispersion dia- 
gram is illustrated in Fig. [6j 

These two dispersion diagrams show a common struc- 
ture in the fluctuations in the form of a band signature of 
the intensity near the gyro-frequency and at higher har- 
monics in the magnetic and compressive dispersion. We 
believe that this is a typical indication for the occurrence 
of ion-Bernstein waves^^M A linear branch with weak 
power is observed at u/k = 1 in the perpendicular disper- 
sion plots of the magnetic field and density fluctuations. 
It corresponds to the linear fast-mode wave in a low-beta 
plasma. A linear Alfven wave does not propagate perpen- 
dicular to the background magnetic field, and therefore 
can be excluded from the interpretation of this branch. 
A coupling between it and the ion-Bernstein modes is not 
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FIG. 5. Dispersion relation of magnetic field fluctuations per- 
pendicular to the background magnetic field at t « 500 for 
Run A. 
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FIG. 6. Dispersion relation of density fluctuations perpendic- 
ular to the background magnetic field at t « 500 for Run A. 
The fluctuations coincide well with the magnetic fluctuations. 



observed, and may presumably not be detectable given 
the spectral resolution of the numerical analysis applied. 

There is further power distributed in compressive and 
magnetic structures at u « 0. These signatures cannot 
be explained as ion-Bernstein waves, because they do not 
have the minimum frequency of Q p and are merely spatial 
structures constant over time. To analyze their nature, 
the possible correlation between the magnetic pressure 
Pb = SB 2 /(8tt) and the density fluctuations Sn can be 
applied. It is defined as 



C 



(SnS\B\ 



(Sn 2 )(S\B\ 



(7) 



where the brackets indicate a certain way of averaging. 
This averaging is here done only perpendicular to the 
magnetic field and in the time domain. Averaging over 
a long time scale corresponds to structures with low fre- 
quency and averaging over short time scales corresponds 
to higher frequencies. The same is true in the spatial 
domain. 

The result of this calculation is shown in Fig. [7) At 



S 5 10 I 

CO 
CD 

co 8 - 




4 6 8 10 

Spatial averaging 

FIG. 7. Correlation coefficient C between fluctuations of 
magnetic pressure and density depending on the averaging 
in space and time for Run A. A strong anti-correlation is 
found up to long averaging times, which correspond to low 
frequencies in the dispersion diagrams. This is an indication 
for the existence of PBSs. The units are given in the above 
introduced normalization. 



high values for the time averaging and at low values for 
the spatial averaging, a strong anti-correlation between 
SPb and Sn is found. This indicates the existence of 
pressure-balanced structures (PBSs), which correspond 
to steepened-up slow-mode waves in the perpendicular 
direction^. They are characterized by their direction of 
propagation and the strong anti-correlation between SPb 
and Sn. The classical slow- mode wave does not propa- 
gate at 90 degrees, but in this direction transforms into 
a tangential discontinuity, which some of the numerical 
structures may represent. Therefore, it seems useful to 
analyze the perpendicular correlation only, since a paral- 
lel slow-mode component is not expected to survive but 
to undergo strong Landau damping. However, it is also 
important to state that our simulation results cannot be 
unequivocal on this issue, because at other positions and 
different time intervals the anti-correlation is not always 
pronounced, and some cases even show a positive corre- 
lation. These findings should be understood just as an 
indication for the presence of PBSs. In this special case, 
the correlation remains positive for only a few averaging 
steps in time. 

The temperature of the treated protons stays mainly 
constant in this case. Also a potential temperature 
anisotropy is not observed. 



B. Results for simulation Run B 

Run B with a higher initial amplitude is in general 
consistent with the previously discussed Run A. There- 
fore, we limit our description to the major similarities 
and differences between the two runs. The magnetic field 
spectrum is shown in Fig.[8j The cascade is also preferen- 
tially oriented in the perpendicular direction. However, 
a larger wavenumber range is filled isotropically around 
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FIG. 8. Two-dimensional power spectral density of magnetic 
field fluctuations in arbitrary units for Run B. The back- 
ground magnetic field is oriented along the k z -axis. Also in 
this case, the cascade is preferentially perpendicular. 
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FIG. 10. Dispersion relation of magnetic field fluctuations 
parallel to the background magnetic field at t ~ 500 for Run 
B. The overplotted red lines indicate again the theoretical 
cold-plasma dispersion relations. 
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FIG. 9. Power spectral density of fluctuations in the magnetic 
field and density at t — 500 along the perpendicular direction 
for Run B. 



k = 0. The compressive component is also very similar 
to the one in Run A (not shown here). 

A cut along the perpendicular direction in the power 
spectra for the magnetic field and density fluctuations is 
shown in Fig. [9j In this case, the spectrum follows more 
closely a slope with the power index of —5/3. 

The dispersion analysis of the parallel fluctuations is 
shown in Fig. [To] and reveals also the same normal mode 
structure as in Run A, yet with a stronger intensity. 

A substantial difference occurs in the perpendicular 
dispersion, which is shown in Fig. 11 The ion-Bernstein 
signatures, 



which have been observed in Run A, are not 
any longer visible in this dispersion diagram and are also 
not present in the compressive dispersion (not shown 
here). Instead, the perpendicular fast mode is more pro- 
nounced. 

The low-frequency part also shows an anti-correlation 
between SPb and Sn under certain conditions and, there- 
fore, seems to represent the same structures as observed 
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FIG. 11. Dispersion relation of magnetic field fluctuations 
perpendicular to the background magnetic field at t « 500 
for Run B. 



in Run A. 

While the wave power in the low-amplitude case was 
too low to produce significant proton heating, this sit- 
uation changed in the case of Run B. The temperature 



evolution of the protons is shown in Fig. 12 Their tem- 
peratures increase, and a temperature anisotropy with 
higher perpendicular than parallel temperature with re- 
spect to the background magnetic field is found. The 
perpendicular temperature raises faster than the parallel 
one and reaches a saturation value of about 1.7 of the 
initial perpendicular temperature. 



IV. DISCUSSION AND CONCLUSIONS 

According to our numerical simulations, the turbu- 
lent energy cascades preferentially into the direction per- 
pendicular to the background magnetic field, which is 
consistent with the recent results of other numerical 
simulations^ 3 -*^ and space observation d 12 * 36 1 3 -^. However, 
there is also a parallel cascade present, and the nature 
of all the relevant fluctuations shall now be discussed in 
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FIG. 12. Temperature evolution of the protons for Run B. 
The temperatures versus time in the parallel and perpendic- 
ular direction with respect to the background magnetic field 
are shown. They are normalized to their initial values. 



dependence upon the direction of propagation. 

The parallel fluctuations seem to be well described as 
a superposition of normal modes. In the range above 
the ion-cyclotron scales they are mainly F/W w aves, as 
it has been suggested before by many authors^- 9138 * 391 . 
This finding is in agreement with observations in the so- 
lar wind indicating that left- and right-handed normal 
modes coexist until a wavelength of the order of l/^ p , 
where the left-handed waves undergo resonant wave- 
particle interactions^] with cyclotron absorption. At 
higher wavenumbers, only right-handed waves can sur- 
vive the transition into the intermediate dissipative range 
of solar wind turbulence^. Linear wave damping seems 
to dominate the dissipation compared to nonlinear damp- 
ing effects, as it has been discussed before in the context 
of interstellar medium heating 40 . 

Parallel A/IC waves are the most prominent of the 
left-handed waves that are known to undergo strong ion- 
cyclotron damping under certain conditions^. The ob- 
served perpendicular solar wind heating can, thus, be 
largely explained by absorption of this left-handed nor- 
mal mode, which is mostly observed in parallel propaga- 
tion. This is most probably also the explanation for the 
establishment of the observed temperature anisotropy 
in the simulation Run B with the higher initial ampli- 
tude. Quasi-perpendicular ion-Bernstein waves may also 
be able to heat the plasma due to cyclotron-resonance 
effects. Therefore, also the here identified ion-Bernstein 
modes could provide a heat source for the ions. The effi- 
ciency of the dissipation of ion-Bernstein waves strongly 
depends on beta and is higher for larger beta values. This 
effect is also discussed by Markovskii et al! 3 -^ in the con- 
text of a cascade of F/W waves. This may be a possible 
explanation for the vanishing of the ion-Bernstein modes 
in Run B because the temperature increase due to the 
resonant heating corresponds to a change in the plasma 
beta. 



The reason for the flatter power index in Fig. [3] com- 
pared to the Kolmogorov scaling, which is instead found 
in Run B as seen in Fig. |9j may be that in the Run A 
with lower initial stirring amplitude the cascade develops 
more slowly, and the turbulence is not yet fully evolved 
at the applied integration time. 

Fast waves themselves can of course also be dissi- 
pated by ions if they propagate obliquely^. This effect 
may, however, be slow when compared to the cyclotron- 
resonant absorption of A/IC waves and needs to be accu- 
mulated over a longer solar wind travel time to become 
significant. The intensity of the ion-Bernstein bands is 
more pronounced for higher electron betas, which is an 
indication for the electrostatic character of these wave 
structures. A kinetic micro-instability, which is able to 
excite ion-Bernstein waves in a way consistent with the 
wave structures observed in magnetospheres, has recently 
been treated in detail with particle-in-cell simulations 4 ^. 

The nature of the perpendicular low-frequency fluc- 
tuations cannot be uniquely identified. There is evi- 
dence for the existence of pressure-balanced structures 
(PBSs), which show the typical anti-correlation between 
5Pb and 5n. It is observed in the solar wind that this 
anti-correlation dominates on shorter time averaging. 
A positive correlation dominates on longer time-scales, 
which is interpreted as the indication for co-rotating in- 
teraction regions as a result of interactions between differ- 
ent solar wind streams with high and low outflow speeds. 
The correlation shows a typical spatial dependence on 
large scales in the solar wind. A positive correlation is 
built up inside 0.7 to 0.8 AU, whereas the anti-correlation 
is already observed closer to the Sun prevailing over a 
large distance ranged Cluster observations also show 
PBSs on smaller scales than the typical low- frequency 
MHD range 4 ^. The origin of these structures, however, 
is unclear. Part of these structures may be generated 
by a nonlinear cascade of the turbulence into the inter- 
mediate wavenumber range, as it is revealed by the sim- 
ulations. Other possible non-propagating perpendicular 
wave structures are mirror modes or Weibel modes. 

The predominance of an anti-correlation in 5Pb and 
5n at some positions supports the interpretation of these 
structures in terms of PBSs. But they are of dynamic 
nature and a minor component driven by the perpen- 
dicular incompressible fluctuations which are the major 
component and represent the key nonlinearly interacting 
degrees of freedom. Therefore, incompressible turbulence 
prevails, but additionally the compressive PBS-like cor- 
relations occur at the same time. The interpretation of 
such fluctuations in terms of static PBSs would be closer 
to a normal-mode picture of the modeled fluctuations but 
perhaps not quite adequate. However, the present simu- 
lation results do not permit a definite conclusion on the 
nature of the weakly compressive fluctuations. 

In a recent model 45 , the transition from strong to weak 
turbulence at small scales is discussed. The authors as- 
sume a completely suppressed parallel cascade for the 
kinetic Alfven waves considered in their model, deduced 



8 



from the suppressed parallel cascade of weak turbulence 
in the incompressible MHD limit. According to their 
model, the anisotropy in the strongly turbulent domain 
is expected to follow the phenomenology of critically bal- 
anced turbulence, which is not applicable to weak tur- 
bulence. Our simulation shows, however, that the paral- 
lel cascade is not fully suppressed, even not in the limit 
of weak turbulence. Yet, the spectral transfer is clearly 
anisotropic and shows a preference for the perpendicular 
cascade. 

The MHD modes with low wavenumbers mainly keep 
their amplitude level during their temporal evolution and 
stay mostly isotropic. A longer integration time might 
also show a preferred direction of the cascade at lower 
wavenumbers, as it was observed by other authors^. 
Waves with higher wavenumbers, however, may hand 
over their energy more easily. This underlines that 
the spectral transfer occurs also in this case locally in 
wavenumber space, as it seems typical for a turbulent 
cascade^. This effect is also incorporated in models con- 
cerning the coronal heating problem, where a nonlinear 
local cascade is often effectively described a s an a dvection 
and diffusion process in wavenumber space^^. 

Normal modes can also be directly excited from the 
ubiquitous thermal fluctuations in a plasma 50 . However, 
their amplitudes usually stay on a very low level, and 
they do not increase with time at the expense of energy 
drawn from lower wavenumbers, as it is observed in our 
simulations. Therefore, it is reasonable to interpret the 
observed wave structures at higher frequencies as being 
the products of processes driven from the low-frequency 
side in the sense of a nonlinear mechanism, instead of be- 
ing of purely thermal origin with a typical energy content 
of the order k^T. 

If the wave power in the inertial range is much higher 
than assumed here, other nonlinear couplings might play 
a role also in the kinetic regime. This can possibly de- 
stroy the normal mode superposition even above k = 1, 
and then lead to a completely different picture with re- 
spect to both the spectral transfer to higher wavenumbers 
and the dispersion structure in the dispersive spectral 
range. 

Parashar et alJ^have shown that the properties of the 
turbulence driver are very important for the generation 
and the further evolution of plasma turbulence. The de- 
pendence on the driving frequency is crucial in the con- 
text of AC coronal heating models, since footpoint mo- 
tions on different timescales are expected to be the main 
source for the turbulence, which is then assumed to heat 
the coronal plasma in these models. 

Finally, our simulation box is small compared to all 
global structures in the solar wind, and the numerical 
conditions are, hence, homogeneous in the above consid- 
erations. There are, however, recent indications that in- 
homogeneities foster the thermalization of wave energy 5 -^. 
This effect may play also a role for the global evolution 
of the solar wind. 
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